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Abstract 

The total cross section of the reaction pp — > ppK^K~ has been measured at excess 
energies Q = lOMeV and 28MeV with the magnetic spectrometer COSY-11. The 
new data show a significant enhancement of the total cross section compared to 
pure phase space expectations or calculations within a one boson exchange model. 
In addition, we present invariant mass spectra of two particle subsystems. While 
the K^K^ system is rather constant for different invariant masses, there is an 
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enhancement in the pK system towards lower masses which could at least be 
partially connected to the influence of the A(1405) resonance.. 
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1 Introduction 

The strength of the kaon-antikaon interaction appears to be very essential 
with respect to different physics topics. It is an important parameter in the 
ongoing discussion on the nature of the scalar resonances oq and /o in the 
mass range of ~1 GeV/c^. Besides the interpretation as a qq meson [1], these 
resonances were also proposed to be qqqq states [2], KK molecules [3,4], hybrid 
gg/meson-meson systems [5] or even quark- less gluonic hadrons [6]. Especially 
for the formation of a molecule, the strength of the KK interaction is a crucial 
quantity and it can be probed in the KK production close to threshold [7] . 
Due to the unavailability of kaon targets for the analysis of KK scattering, the 
kaon pair production in multi particle exit channels like pp ppK^K~ is the 
only possibility to study this interaction by selecting the appropriate kinematic 
region of the phase space distribution. Besides the KK subsystem, informa- 
tion about the KN system is of equal importance especially in view of the 
actual discussion on the structure of the excited hyperon A(1405) which is 
considered as a 3 quark system or a KN molecular state [8]. Up to now the 
scattering length a^-p has been mainly determined on kaonic hydrogen. But 
the situation is not yet clarified since first, the results of former measurements 
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[9,10,11,12] and preliminary results at DEAR [13,14] are in disagreement and 
second, it has been shown that contrary to pionic hydrogen, the isospin vio- 
lating correction cannot be neglected in the kaonic case [15]. Due to these still 
open questions and the fact that the analysis of former KN data (cf. figure 22 
in [12] and references therein) have a rather large variation, new low energy 
pK' scattering data can provide an independent contribution to this impor- 
tant issue. 

Furthermore, a precise knowledge of the cross sections and a good under- 
standing of the kaon and antikaon interaction with the nucleon is an essential 
ingredient for calculations of medium effects [16] being related to open ques- 
tions of astrophysics [17]. This is because in dense matter processes there are 
secondary production mechanisms on hyperons such as t:Y — > K'N making 
it necessary to understand the production above threshold [18]. 
While the database on the elementary creation covers a wide energy range 
[19,20,21,22,23,24], low energy data on the K~ production are less available 
[25,26,27]. In the near threshold regime, the excitation function might show a 
significant difference compared to the expectation of a pure phase space be- 
cause final state interaction effects are predominant at low relative energies of 
the outgoing particles [28]. 

Due to the mentioned aspects together with the tendency becoming appar- 
ent that the available data seem to lie above theoretical expectations, we 
performed two new measurements of the total cross section of the reaction 
PP ~^ ppK^K~ at excess energies of Q = lOMeV and 28MeV [29] in order to 
further study this enhancement and its strength. In the next section, we de- 
scribe the experimental technique followed by the presentation of the results. 



3 



2 Experiment 



The measurements of the pp — »■ ppK^K~ reaction were performed with the 
internal experiment COSY-11 [30] at the COoIer SYnchrotron COSY [31] in 
Jiihch with a beam momenta oip = 3.333 GeV/c and 3.390 GeV/c correspond- 
ing to excess energies of Q = 10 MeV and 28 MeV, respectively. Both energies 
are below the 0(1020) meson production threshold. The detector which is 
shown in figure 1 is designed as a magnetic spectrometer. Using a hydrogen 



Fig. 1. The experiment COSY-11 with the main components. The overlayed 
box shows a schematic view for an exemplary event of the reaction channel 
pp — > ppK^K^. The sudden stop of the kaon track in the three dimensional picture 
indicates its decay. 

cluster target [32] in front of one of the COSY dipole magnets, positively 
charged particles in the exit channel are bent more - compared to the remain- 
ing beam protons - towards the interior of the ring where they are detected in 
a set of three drift chambers [33]. Tracing back the reconstructed trajectories 
through the known magnetic field [34] to the interaction point allows for a 




scintillator S1 
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momentum determination. In combination with a subsequent time of flight 
measurement over a distance of 9.4 m between two scintillation hodoscopes SI 
and S3, these particles are identified via their invariant masses. Due to the 
decay of the kaon, the probability that it reaches the stop counter S3 is in 
the order of a few percent. Therefore, an indirect reconstruction of the time 
of flight is used. After the determination of the two protons' four momenta in 
combination with the known length of their flight path from the target to the 
SI detector, the time of the interaction is calculated. This time is then used 
as the start for the kaon's time of flight between the point of interaction and 
the crossing in the SI scintillator in order to derive the four momentum of the 
K+. 

Two additional detector components are mounted in front of the dipole mag- 
net close to the target and inside the dipole gap both consisting of a scintillator 
and silicon pads. While the first is used to measure the coincident proton of 
the pp-elastic scattering, the array in the dipole gap serves to detect the K~. 
The analysis for the reaction pp ppK^K^ proceeds in several steps. First, 
events with less than three reconstructed tracks are rejected. For the remain- 
ing data, figure 2 shows the squared invariant masses for those two tracks 
that could be assigned to a hit in the S3 scintillator. A clear separable peak 
for two protons is visible. With the described indirect method for the time 
of flight, the four momentum of the third positive particle is deduced. Figure 
3(a) shows the squared mass of the third particle X+ versus the missing mass 
of the three particle system assuming that X+ is a kaon. The z-axis is in a 
logarithmic scale. Besides the huge horizontal pionic band there is a less pro- 
nounced band structure along the kaon mass. At the cross point of both green 
dashed lines (corresponding to the literature value of the charged kaon mass 
[35]) a separated group of events is visible. 
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Fig. 2. Identification of the two protons via the squared invariant mass for events 
with three reconstructed tracks. Here, the two particles are encountered which reach 
the final scintillator S3. 




missing mass"^ PpK"*" [GeV^/c^] missing mass ppK"*" [GeN^lc^} 

(a) (b) 
Fig. 3. a) Squared invariant mass of the third positive particle X versus the missing 

mass assuming m-^^ = m^+. The green dashed lines indicate the literature kaon 

mass. The z-axis is in logarithmic scale, b) Missing mass of the ppK~^ system for 

events in between the blue solid lines in the left plot. An additional request for 

hits of certain modules in the SI scintillator was required. The shaded spectrum 

includes further cuts described in the text. It is worth noting that the missing mass 

resolution is less than the size of the bins. 

The projection within a broad band around the kaon mass (horizontal blue 
solid lines in fig. 3(a)) is shown in figure 3(b). Here, additional cuts on the 
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segments in the SI scintillator were applied since the protons are passing the 
detector closer to the beam pipe than the K^. These cuts were adjusted with 
Monte Carlo simulations for both energies separately. The missing mass spec- 
trum (fig. 3(b)) shows a clear peak at the K~ mass and a broad physical 
background towards lower missing masses. The latter is understood mainly 
with the excited hyperon production pp — > pi^+A(1405)/E(1385) where the 
second proton stems from the decay of the resonance. Additionally, pion pro- 
duction channels contribute where pions are misidentified as kaons [26]. 
A fit of the kaon peak with a Gaussian and a polynomial function describing 
the background results in a missing mass resolution of ct^m ~ IMeV/c^ for 
both energies. Additional cuts using the detector mounted inside the dipole 
gap can be used to drastically reduce the background. The explicit procedure 
can be found in [26,29] and was performed to cross-check the signal result- 
ing from the K~ production. The principle procedure is to calculate with the 
known K~ four momentum its expected hit position in the silicon pads inside 
the dipole gap and to compare it with the experimental location. While for 
real K~ events this should be the same, a background event does not show a 
correlation of these two values. The final and nearly background free missing 
mass spectrum is shown in figure 3(b) (shaded area). The loss of kaons due to 
its decay agrees perfectly with Monte Carlo studies. Nonetheless, the deduc- 
tion of the total cross section is based on the non-shaded missing mass plotted 
in figure 3(b) due to the higher statistics for the kaon signal. The background 
contribution was extracted assuming a polynomial function to describe its 
shape and the signal to background ratio was 5.4 and 3.8 for Q = 10 MeV and 
28MeV, respectively. 
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3 Results 



3.1 Luminosity and efficiency 

The absolute normalization of the counting rate requires both the knowledge 
of the luminosity £ and the total detection efficiency E including the geo- 
metrical acceptance of the detector and the reconstruction efficiency. For the 
determination of the luminosity, the elastic proton proton scattering is used. 
While one proton is registered in the main detector as presented above for the 
reaction products and therefore its four momentum is determined, the second 
proton is registered in the scintillator and silicon pads in front of the dipole 
magnet [36] . The differential counting rates are normahzed to the EDDA data 
[37]. The determined integrated luminosity / Cdtis given in table 1, including 
statistical and systematical errors. 

The total detection efficiency for the K^K^ production reaction was investi- 
gated using Monte Carlo simulations based on the GEANT 3 code [38]. This 
software package has been designed to completely describe the response of the 
detector. It is important to mention that the COSY- 11 results obtained so 
far at similar excess energies but in other reaction channels are in very good 
agreement with measurements at other laboratories. In particular, results of 
measurements at COSY-11 for the pp pprj' reaction [39] with a beam mo- 
mentum by only 1% lower to the one reported in the present article, are in 
excellent agreement with the results obtained at the SATURNE facihty [40]. 
Using this software package, for each generated event a detection system re- 
sponse is calculated and the simulated data sample is analyzed with the same 
program which is used for the analysis of the experimental data. The to- 
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tal efficiency E for tfie free reaction pp ppK'^K~ including tfie pp final 
state interaction is listed in table 1. Here, the identification of two protons 
and a in selected segments of the SI is required. The systematical error 
comprises the detection and reconstruction efficiencies, the decay of the kaon 
and the variance of the efficiency resulting from the uncertainty of the beam 
(Ap/p < 10~^) that has been extracted to be in the order of 1-2 MeV/c [29]. 
Due to the scaling of the efficiency with l/Q, this translates to a variation of 
8 by 9% and 1% at Q=10 MeV and Q=28MeV, respectively 
Furthermore, we included an estimate of the influence of higher partial wave 
contributions. Prom measurements at even higher excess energies [26] and the 
extracted angular spectra [29] there is no indication of a substantial contri- 
bution from higher partial waves. In consequence, the additional term for the 
systematical error of a few percent is rather overestimating and a conservative 
upper estimate. 

The inclusion of the pp-FSI to the Monte-Carlo studies results in a relative 
change of S by 10% [41,39]. Since the various descriptions for the pp-FSl differ 
from each other by only around 30%, the contribution to the relative system- 
atic error of £ from the inaccuracy of the knowledge of the FSI is 3% (30% out 
of 10%). Possible influences stemming from the pK and KK FSI are negligible 
compared to the pp-FSI ^ and therefore were not included. 



^ For the pp-FSI the scattering length app = 7.8 fm [42] is much larger compared 
to the one for the pK system (a^x) being less than 1 fm. Since the FSI is in first 
order proportional to the squared scattering length this neglect is reasonable. 
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3.2 Total cross section 



Using the knowledge of the integrated luminosity and the overall efficiency, the 
number of events registered for both excess energies can be transformed into a 
total cross section atot- For both energies, we obtained the final results given in 
table 2, where the contributing systematical errors were added quadratically. 
The new results together with previous measurements are compiled in figure 
4(a) together with some theoretical expectations. 

It is obvious that at low excess energies the data points lie significantly above 
the expectations indicated by the different lines that are all normalized to the 
DISTO point ai Q — 114 MeV. The reason for this choice is the fact that 
the effect of the nucleon-nucleon FSI diminishes with increasing excess energy 
since it significantly infiuences only that part of the phase space at which 
nucleons have small relative momenta. While this fraction stays constant, the 
full phase space volume Vps grows rapidly: A change from Q = 1 MeV to 
Q — 10 MeV corresponds to a growth of Vps by more than three orders of 
magnitude. As a result the S-wave pp-FSI is of less importance for higher 
excess energies where it affects a small fraction of the available phase space 
volume only [43] . Additionally, a possible contribution of higher partial waves 

Table 1 

The determined total integrated luminosity including statistical and systematical 

errors and the total detection efficiency for both Q values. 

g = 10MeV Q = 28MeV 

JCdt [pb-^] 2.770 ± 0.045 ± 0.011 2.270 ± 0.064 ± 0.006 
£ [%] 1.238 ± 0.129 0.308 ± 0.027 
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Table 2 

Total cross section for both Q values including statistical and systematical errors. 

Q [MeV] <Jt,t [nb] 

10 0.787 ± 0.178 ± 0.082 

28 4.285 ± 0.977 ± 0.374 

aX Q = 114 MeV should even increase the total cross section at this energy. 
Therefore, their inclusion in the calculations would result in an even stronger 
discrepancy between the calculations and the data at low Q values. 
The pure non-relativistic phase space (dashed line in fig. 4(a)) differs from 
the experimental data by two orders of magnitude at Q = 10 MeV and a 
factor of five to ten aX Q = 28 MeV. In comparison to that, the inclusion of 
the j9p-FSI (red solid line) by folding its parameterization known from the 




1 10 10^ 4 10 10^ 



Q [MeV] Q [MeV] 

(a) (b) 
Fig. 4. a) Total cross section as a function of the excess energy Q for the reac- 
tion pp — > ppK^K^. The former experimental results are taken from [25,26,27] and 
shown together with the new results (triangles) for which data the statistical error 
(horizontal right markers) and the statistical plus systematic errors (horizontal left 
markers) are shown, b) Ratio of the experimental data aexp over the calculation for 
pure phase space and for the parametrization of the pp-FSI. 
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three body final state with the four body phase space is clearly closer to 
the experimental results but does not fully account for the difference ^ . The 
solid line representing the calculation within a one-boson exchange model 
[44] reveals a similar discrepancy as the pp-FSI parameterization. This model 
includes an energy dependent scattering amplitude derived from the fit of the 
total cross sections in K'^p — > K'^p [45] while the pp-FSl was not included, 
yet. Up to now, there is no full calculation available but the new data demand 
further theoretical efforts in order to give a complete picture of the K^K~ 
production. 

The enhancement of the total cross section at low energies could be partly 
induced by the opening of the neutral (K^K^) kaon pair channel with a mass 
splitting compared to the charged {K'^K~) kaon pairs of about 8MeV. Since 
calculations show a substantial influence of the opening neutral channel on 
the TTTT — K^K~ cross section (cf. figure 2 in reference [7]), an effect on the 
excitation function of the reaction pp ppK^K" could be expected at low 
energies as well. A quantitative estimation of this effect has to be calculated 
within a complete coupled channel model taking into account this K^K~ ^ 
K^K^ transition. However, the kinematical situation of the four body final 
state is expected to strongly suppress the effect [46] compared to the two 
body final state tttt — > K'^K~ . 



^ The parametrization presented differs from that in figure 6 of reference [26] . The 
reason is that the formula used in [26] was originaUy derived for the three body 
final state. Our new ansatz [29] avoids this approximation by splitting the integral 
over the four body phase space into two parts whereas one of them contains a three 
body subsystem for which the pp-FSI parametrization is known. 
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3.3 Differential spectra 



The knowledge of the four momenta allows to investigate differential observ- 
ables. Here, in particular the invariant mass of several subsystems can be 
exploited. For two particles i and j in the exit channel, the square of the 
invariant mass^ m^j is given by — (Pj + Pj)^ with Pj being the four mo- 
mentum of the particle i. The events for the following plots were chosen under 
the condition, that the missing mass in figure 3(b) is within the region 
of the kaon (0.235 GeVVc^ < m\ < 0.25GeVVc^)- Furthermore, the events 
were corrected by the detection efficiency and then normalised to the phase 
space distribution . Generally for each mass bin a separate background sub- 
traction should be performed in order to make a reliable interpretation. For 
high statistics this in fact has been done (e.g. in case of the pp pprj [47]). 
For the present data, however, the background contribution is that small (as 
seen in figure 3(b)), that such a procedure could be avoided. 
For a pure phase space distribution, the ratio should be flat as it 

rather is in case of the K'^K' system shown in figure 5(a) for both Q values. 
In a former publication [26] it has been shown, that a possible influence of the 
scalar resonances ao(980)//o(980) is not distinguishable from a pure s-wave 
distribution within the current statistics. 



^ In case of fn^j^j^ or m^^_ , for each event both protons are taken into account so 

that each invariant mass m^j^ has two entries per event. 

^ Therefore, instead of the mass spectrum rriij we show the ratio R^^ := 

where is the number of events registered in a certain mass bin of the subsystem 

ij and /'^ps is the number of accepted Monte Carlo events in the same mass bin 

ij 

generated with an underlying phase space distribution. Deviations from the phase 
space directly reflect into a non flat distribution of Rij. 
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Figure 5(b) shows the double ratio it! := RpK-l ^pK+- 'The normahsation of the 
pK~ system to the pK^ has the advantage, that any systematic errors from 
misunderstood inefficiency should in principle affect both systems in a similar 
way and therefore mainly cancel out. In consequence, since the interaction in 
the pK^ system is known to be rather weak ^ , any non flat distribution in 
this ratio R might better indicate an interaction between the proton and the 
negative kaon than the pure spectrum of m^j^- alone. There is a clear increase 
towards lower invariant masses in the double ratio R. This could result from 
the final state interaction of the two particles or partially be a refiection of 
the A(1405) ^ or a mixture of both. 

The shown distributions might trigger further investigations with improved 
statistics with the new installation of the WASA detector [48,49] which com- 
bines a large acceptance of nearly 47r with the simultaneous detection of 
charged and neutral particles. 

4 Summciry 

The COSY-11 collaboration has extended its studies on the elementary K' 
production by measuring the reaction pp — > ppK^K~ at excess energies of 
Q = lOMeV and 28MeV resulting in total cross sections exceeding the ex- 
pectations for a pure phase space drastically. The measurement is based on a 
^ The experimental distribution for this subsystem is indeed within the error bars 
not deviating from the pure phase space (similar hke in the case of the K^K~ 
system in figure 5(a)). 

^ The S(1385) resonance certainly will influence this system as well. However, due 
to it's lower mass and slightly smaller width [35], it's contribution should be less 
significant than that of the A(1405). 
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(a) (b) 
Fig. 5. a) Invariant mass m^+^_ for both excess energies normalised to the accepted 

MC events generated with a phase space distribution, b) Invariant mass for the 

system pK^ divided by that for the pK^ system. For a better clearness the data 
points at (5 = 10 MeV were slightly shifted right for both pictures. 

kinematically complete reconstruction of the positively charged ejectiles while 
the negative kaon is identified via the missing mass. An absolute normaliza- 
tion of the counting rate is achieved via a simultaneous measurement of the 

pp clastic scattering. 

The new results for the total cross section in the reaction pp — > ppK^K~ are 
clearly showing that towards the lower Q values the data are exceeding any 
expectations both from pure phase space with and without the pp FSI and 
a calculation within a meson exchange model. To further study the strength 
of the enhancement of the total cross section at low Q values, the COSY- 11 
collaboration will remeasure the data point at Q = 6 MeV [50] in order to 
significantly reduce the statistical error. 

Within the limited statistics the differential m^j^- distribution normalised to 
the pK^ system shows a pK~ interaction which might have a connection to 
the A (1405). It is too early to extract quantitative couphng strengths infor- 
mation from the present data on the hyperon resonances S(1385) and A(1405) 
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where especially the structure of the latter one is under discussion. But the 
present data clearly demonstrate the sensitivity in this ppK~^K~ final state 
to the KN interaction which is an important ingredient in the interaction of 
the A(1405) as a bound KN system. 
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